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Abstract 

Our analysis shows how the covariant chaotic behavior characterizing 
q-i the evolution of the Mixmaster cosmology near the initial singularity can 

5^ | be taken as the semiclassical limit in the canonical quantization performed 

by the corresponding Hamiltonian representation. 

> 

X 

1 Cosmological Framework 

Though the Standard Cosmological Model (SCM) is based on the highly sym- 
metric Friedmann-Lemaitre-Robertson- Walker (FLRW) solution of the Einstein 
equations, nevertheless don't exist neither theoretically neither experimentally 
any evidences to prevent that our universe underwent a more general dynamics 
in the very early stages of its evolution, and only in a later phase it isotropized 
reaching a complete agreement with actual experimental data. 
The simplest generalization of the FLRW dynamics consists of the so-called 
Bianchi models, whose anisotropic evolution, for some of them, can be repre- 
sented as a FLRW model plus a gravitational waves packet fl], Among this 
classification, the types VIIlQ and IX appear as the most general ones allowed by 

1 All the considerations we will develop for the type IX apply also to the VIII one since, close 
to the singularity, they have the same morphology. 
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the homogeneity constraint^ and their asymptotic evolution toward the cosmolog- 
ical singularity manifests a chaotic-like behavior @. The cosmological interest 
lies in the IX model (the so-called Mixmaster Q), which has the same space 
symmetries as the closed FLRW universe, and whose dynamics allows the line 
element to be decomposed as 

ds 2 = ds 2 - 5 {am Git m dx t dx k (1) 

where dso denotes the line element of an isotropic universe having constant pos- 
itive curvature, is a set of spatial tensors^] and #( a ) (&)(£) are amplitude 
functions, resulting small sufficiently far from the singularity. 
We dedicate our analysis to find a precise relation between the Mixmaster deter- 
ministic chaos and the quantum behavior characterizing the Planckian era, show- 
ing how the invariant measure for the former, provided general Misner-Chitre-like 
coordinates (MCI) J|, ||, is independent of the time gauge and coincides with 
the stationary probability for the semiclassical limit of the latter. 



2 Billiard Representation 

Using generic MCI variables (£, 9, r) the dynamics is described by the two- 
dimensional canonical variational principle 

S J (M' + Pe9' - FHadm) dv = 0, (2) 
where / is a generic function, 

H ADM = V^+U, e 2 = gV + T ( 3 ) 

q z 



where [/(£, 9, rf) denotes the corresponding potential term and q = g(£) = \/£ 2 — 1. 
Moreover the equation for the temporal gauge reads 

12D o f df , 
N{rj) = - e^j-r' (4) 

n-ADM CLT 

so that our analysis remains fully independent of the choice of the time variable 
until the form of / and r' is not fixed. 



2 IX's geometry is invariant under the SO{2>) group. 
3 These tensors satisfy the equations 



Q(a)(b);l _ _/ n 2 _ 3")£<( a )( b ) g( a )( b ) fe _ q g( a )( b )'» _ q 

in which the Laplacian is referred to the geometry of the sphere of unit radius. 
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For the following developments it is of key interest the following relation 

d (Hadm/') 9 {HadmI') 



(5) 



drj dr] 

The function / (rj) plays the role of a parametric function of time and actually the 
anisotropy parameters 00 H (i = 1,2,3) are functions of the variables (£, 9) 
only 0. In the domain T# where all the Hi are simultaneously greater than 0, 
the potential term U can be modeled by the potential wallsQ 

= Goo {H x ) + e M (H 2 ) + e M (H 3 ) (6) 

therefore, by (^), in the ADM Hamiltonian becomes (asymptotically) an 
integral of motion 

f &Hadm dE = 

v^.nerJ df " df (7) 

( Hadm = V^+U = e = E 
where E is a constant. In view of (|?]) the variational principle (§) reduces to 

^ / (Pfdf + Pe^ - #df) = ^ / (p«^ + P*d0) = . (8) 

Following the standard Jacobi procedure to reduce this principle to a geodesic 
one we get, for the closed domain region r#, the Riemannian line element 



ds 2 = E 2 



de 



, 2 _ T + (^ 2 -l)^ 2 J. 0) 

Since the above metric (|9|) has curvature scalar R = — the point-universe moves 
over a negatively curved bidimensional space on which the potential wall (|6]) cuts 
the region indeed the invariant Lyapunov exponent for the dynamical flux 
associated to @ reads A„ = 1/E > |. The point-universe, bouncing on the 
potential walls, is reflected from a geodesic to another one, making each of them 
unstable. By itself, the positivity of Lyapunov number is not enough to ensure the 
system chaoticity, since its derivation remains valid for any Bianchi type model, 
but for the Mixmaster case the potential walls reduce the configuration space to 
a compact region (r#), ensuring a real chaotic behavior. 



3 Statistical Mechanics Approach 

For a Statistical Mechanics reformulation of the dynamics, we adopt in (Q) the 
restricted time gauge r' = 1, leading to the variational principle 

5 1 \ n Tf + Vd Tf ~ Hadm ) df = °' (10) 



4C 



6oo (a;) 



+oo if x < 
if x > 
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In spite of this restriction, for any assigned time variable r (i.e. rf) there exists 
a corresponding function / (r) (i.e. a set of MCI variables able to provide the 
scheme presented in Section ||) defined by the (invertible) relation 

dr 12D K ' K J 

Hence the analysis to derive the invariant measure for the system follows the same 
lines presented in || f| . Indeed we got a suitable representation of the Mixmaster 
chaoticity in terms of a two-dimensional point-universe randomizing within r#, 
admitting an "energy-like" constant of motion e = E, then well-described by a 
microcanonical ensemble, whose Liouville invariant measure reads[] 

dg = A5 (E — e) d^d9dp^dp e , A = const.. (12) 

After the natural positions 

Pf = -cos0, po = £qsm(J), (13) 
Q 

being < <p < 2ir, and the integration over all values of e[], we arrive to the 
uniform invariant measure |TD|, |5| H 



dfi = Woo (f , 6, 4>) d^ddd<p = -J-d£d9d0 . (14) 
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The key point of our analysis is that any stationary solution of the Liouville 



theorem, like (14), remains valid for any choice of the time variable r; indeed in 
|| the construction of the Liouville theorem with respect to the variables (£, 9, </>) 
shows the existence of such properties even for the invariant measure (|l~4l). 

More precisely, in agreement with the analysis presented in ||, during a free 
geodesic motion the asymptotic functions £ (/) ,9 (/) and (/) are provided by 
the simple system 

, d6 sin0 d(j> £sin</> 

— = gcos0, — = , — = (15) 

df df q df q 

and therefore over the reduced phase spaceQ {£, 9} ® the distribution 
behaves like the step-function 

«-({,»,« = ( 8^ V Wl^**** (16) 

o v {?,»,<#>} 0r„®si 



5 <5 (x) denotes the Dirac function. 

6 The dependence on the initial conditions doesn't contain any information about the system 
chaoticity. 

7 denotes the ^-circle. 
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stationary solution of the Liouville theorem 



sin09woo ^sin^^Woo 

gcos0— — + — H1T = - 17 

ot, q o9 q d(p 

If now we restrict our attention to the distribution function on the configuration 
space T H 

q(£,6)= rw oa (C,9A)d(f>, (18) 
Jo 

by fllTD we get for such reduced form the two dimensional continuity equation 



qCOS ^ + —^f = ° (19) 

and the microcanonical solution on the whole configuration space {£, 9} then 
reads 

Qoo(U) = i v ^> Gr * . (20) 

{ o v {£,0}^ 

4 Quantum Origin of the Chaos 

The main result of the above Section |3| [0, § is the proof that the chaoticity of the 
Bianchi IX model above outlined is an intrinsic feature of its dynamics and not an 
effect induced by a particular class of references. Appearing this intrinsic chaos 
close to the Big Bang, we infer that it has strict relations with the indeterministic 
quantum dynamics the model performs in the Planckian era. The link between 
quantum and deterministic chaos is searched in the sense of a semiclassical limit 
for the canonical quantization of the model. 

Indeed the asymptotical principle (H) can be quantized by a natural Schroedinger 
approach || 

dip 

^~q~ = KadmiP , (21) 

being if) = if)(r,£,9) the wave function for the point-universe and, implementing 
*Hadm (see (Q)) to an operator^, i.e. 

d d 
Pi -> Pi = ~ iTi Q£ ' Pe -* Pe = ~ ih Q§ ' ( 22 ) 



3 The non vanishing canonical commutation relations are 

£,P£ =ih, 9,pg =ih. 
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the equation fl2"T| ) rewrites explicitly, in the asymptotic limit U — > Uoo, 

where we left Uoo to stress that the potential cannot be neglected on the entire 
configuration space {£, 9} and, being infinity out of Th, it requires as boundary 
condition for if) to vanish outside the potential walls ip (dT H ) = 0. Since the 
potential walls Uoo are time independent, a solution of Q23D can be taken in the 
form 

CO 

V>(T,£,0) = £c n e-^/V(£,0) (24) 

n=l 

where c n are constant coefficients and we assumed a discrete "energy" spectrum 
because the quantum point-universe is restricted in the finite region T H . Once 
taken an appropriate symmetric normal ordering prescription, the squared rela- 
tion from (^) and the position ( P^) lead to the eigenvalue problem 



d d I d I d 
~ q d£ q d£ ~~ qd6qd6 



E 2 -U E 2 

± -'n w co ^oon /ocr\ 

= J2 <Pn = —tf-^n- (25) 

The quantum equation (|23|) is equivalent to the Wheeler-DeWitt one for the same 
Bianchi model, once separated the positive and negative frequency solutions, with 
the advantage that now r is a real time variable. 

In what follows we search the semiclassical solution of this equation regarding 
the eigenvalue Eoo n as a finite constant (i.e. we consider the potential walls as 
finite) and only at the end of the procedure we will take the limit for Uoo ©• 
We infer that in the semiclassical limit when H — > and the occupation number 
n tends to infinity (but rih approaches a finite value) the wave function (p n — >• if 
and Eoon — > ^oo so that we have 




<P (£, 0) = (e, 9) exp z^^- , (26) 



where r and S are functions to be determined. 

Substituting ( p6|) in ( p5|) and separating the real from the complex part we get 
two independent equations, i.e. 



classical term 



0(l/ft) 
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The dominant term in (|27| ) reduces to the Hamilton- Jacob i equation and its 
solution can be easily checked to bef] 



(29) 



where k is an integration constant. Through the identifications 
dS dS f 

~Q£ = P Z ' ~QQ = Pe S = J (Pt d ^ + Pdde "> ' ( 30 ) 

( p7|) is reduced to a mere algebraic constraint which is the asymptotic one TC\ DM = 
E 2 = and is solved by fll3|) replacing e = E^, whose compatibility with ( p9|) 
and ( PO] ) is then obtained using the equations of motion fli"5D which provide 



e 



ctgip =>- J£ 2 — 1 sin ip = c , (31) 



where c = const.; the required compatibility comes from the identification k = 
EooC. Since E^ — > ioo outside Th the solution <p (£, 0) vanishes out of the billiard. 
The substitution in ( p8|) of the positions ( Jl3| ) leads to the new equation 

, <9r sin (j) dr 
gcos0— + ~ = 0, 32 

which coincides with (|19D, provided the identification r = g^; this result ensures 
the correspondence between the statistical and the semiclassical quantum analy- 
sis. 
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